AD-A283  058 


SuperMariah: 

A  Similarity-based  Method  for  Determiniiig  Wind,  ^ 
Temperature,  and  Humidity  Profile  Structure  in  the 
Surface  Boundary  Layer 

by  Henry  Rachele 
Frank  V.  Hansen 
Arnold  Tnnick 

Battlefield  Environment  Directorate 
LisaMangoso 

New  Mexico  State  University 


ARL-TR264  June  1994 


% 


94-24207 


rr? 


Appmued  far  pub^itkase;  distribution  Is  unUmOetL 

94  8  01  Oig 


M0TZCB8 


Diaelaiji«rs 


The  fiadiaga  la  this  raport  mra  aot  to  bo  ooaotxuod  as  aa  offied 
D^artaoat  of  tbo  Araj  pesitioa#  ualooo  so  dooigaatod  by  ett 
autborlaad  doenBoats. 

'Tf 

Thm  oitatioa  of  trada  aaaaa  aad  aaaaa  of  aaattfaotaxars  la  ti 
raport  is  aot  to  ba  ooastruad  aa  offleial  Oovazaaaat  iadorsaaoat: 
approval  of  ooaaareial  pxodnets  or  sarvleas  rafaraaead  baraia»~ 


'fi 


Oaatruetloa  Hotlea  .'-M 

Wbaa  tbia  doouaaat  la  ao  loagar  aaadad^  daatroy  it  by  aay  aaa 
that  will  pravaat  disolosuro  of  its  ooataats  or  raooastraotioi^ 
tba  dooaaaat.  li 


REPORT  DOCUMENTATION  PAGE 

Form  Appnrmd 
(3MB  No.  OTOB^m 

rubbc  nooning  bufOtn  tar  Um  coMcnen  of  mtannotion  n  titmwttO  to  ororogc  <  Mur  ptr  mponio.  inclMHng  tht  tuiw  tor  wwounng  imiructiora.  MorOwig  onduig  data  iouron. 
gotMtino  and  momuimng  IM  dou  noidod.  and  comotaong  and  raviaunng  tha  cotlartion  of  information,  tand  commtno  ragarding  tfm  burdan  auimaia  or  any  otlwr  tipan  of  ttaa 
coHacoon  of  information,  including  tuggattiont  tar  raduong  tlm  burdan.  to  Waafungton  naadguanan  larvicaa.  Oiractorata  tar  information  Oparatiam  and  Naporta.  121$  iotfanon 
Oa«nHighway.$uital204.  ariington.  VA  22202-4302.  and  to  tna  Offica  of  Managamant  and  tudgat.  Oapaiwork  daductioo  froiact  (070a-om>.  dfailungten.  OC  20$03. 

4.  TITU  AND  SUITITLf 

SupariMarith:  A  StmilariQr-baaed  Mathod  for  Datemiiiiim  IX^nd, 

Tanpaialiire,  and  Humidity  Frafik  Stniduie  in  the  Suiftoe 

Boundaiy  Lqrer 

5.  FUNDING  NUMBERS 

«.  AUTHOmS) 

Hamy  Reehela,  Fiank  V.  Hanaan,  and  Anudd  Tkmidc  (ARL) 

Liaa  Manguao  (New  Meadao  State  Univanity) 

7.  PERFORMING  ORGANUATION  NAME(S)  AND  AOORESS<ES) 

U.S.  Amy  Raaeaidi  Labotatny 

Battlefield  Enviiaoment  Oinetanta 

Attn:  AMSRL-BE-W 

White  Sanda  Miaaile  Range,  New  Meacieo  88002-3501 

B.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ARL-TR-264 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS<ES) 

U.S.  Amy  Raaeaiah  Labontoiy 

2800  PoiMier  Mm  Road 

Ada^ihi.  MD  20783-1143 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

ARL-TR-264 

11.  SUPPLEMENTARY  NOTES 

12a.  OISTRWUTION/AVAILARILITY  STATEMENT 

/^piovad  for  public  laleaae;  dkcribiilian  ia  unGmilad. 

12b.  DISTRIBUTION  CODE 

A 

13.  ARSTRACT  (Maximum  200  words) 

Siyeriitoriah  ia  a  new  methodology  gar  aatinintiag  the  auifaeelayaraimilBiiQfacalingwaialarta  for  mind,  tempeiatuia,  and  apecifie  humidily. 
uaing  data  firom  totwar-boma  aenama.  It  ia  conceptually  baaed  on  tha  Dynamic  Siniilari^y«rfPh>$wTlieoiy.  We  diacum  the  origin  and  logic 
of  the  Mariah  mpioach  (tha  piadecaaaor  <rf  SiyetMarinh)  and  then  p recent  tha  cflanaion  of  Mariah  nauhing  in  SupciMaxiah.  ThiaiqMft 
lavietvB  the  con^leta  eat  of  modal  aquathmi  and  piaaanta  aavanl  aaanylm  to  iHuatiata  aompariaona  ot  the  Mariah,  SupeiMariah,  and  the 
more  treditional  O’KEYPS  fennulatiooB. 

14.  SUBJEa  TERMS 

nucremeterology,  ObuUwv  aimilarity,  turfKa  boundaiy  Iqrer, 

O’KEYPS  function,  tcmpeiature,  ^lecific  humidi^,  wU  apaed 

15.  NUMBER  OF  PAGES 

33 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclaaaifiad 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclaaaifiad 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclaaaified 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-01 -280-5500  Standard  Form  298  (Rav.  2-89) 


rmcnbtd  br  *mi  Std  23<-ia 

2m-102 


Contrats 


1.  IntroductkMi .  3 

2.  Mariah .  S 

3.  SuperMuiah . 9 

4.  The  SuperStar  Optioii .  13 

5.  Discussion .  17 

6.  Conclusions . 25 

References . 27 

Bibliography . 29 

Appendix  A.  Obukhov  Length  for  Stable  Conditions . 31 

Distribution . 35 


figures 


1.  Si^rMariah  calculated  profiles  compared  to  data  from  REBAL  ’92: 

9  July  1992  1000  CST  unstable  flow,  (a)  Temperature  where  0*  =  -.4856, 

(b)  Wind  speed  where  u*  =  .544,  (c)  Spedfic  humidity 

where  q*  =  -5.335  *  lOe-5 . 

2.  SuperMariah  calculated  profiles  compared  to  data  from  REBAL  ’92: 

15  July  1992  2245  CST  stable  flow,  (a)  Temperature  where  0*  =  .0742, 

(b)  Wind  speed  where  u*  —  .184,  (c)  Specific  humidty 

where  q*  =  1.474  •  lOe-5 . 


Aoo«sslon  For  j 

BTIS  0RA4I 

DTIC  TAB 

n 

UnaRr.o’oncod 

□ 

J»s  t If : c  at i on — 

By - 

DlsWib«tl?n/ . 


Availsbillty  CodfiS_ 
tA'rrUi  r nd/or 


Diet 


V  2lt- 


11 


12 


1 


3.  Comparison  of  SuperMariah  and  SuperStar  profiles,  N^raska 
7  August  1956  1205  CST  (Lettau,  1957).  (a)  Potential  temperature 
where:  SuperMariah  fT  =  -.3989,  Superstar  ^  =  -.4024, 

(b)  Wind  speed  where:  SuperMariah  u*  =  .4741, 

SuperStar  u*  =  .4578,  (c)  Specific  humidity 

where:  SuperMariah  q*  =  -1.185e-5,  SuperStar  q*  =  -1.117e-5 . 16 

Tables 

1.  Comparison  of  SuperMariah,  Mariah,  and  O’KEYPS  schemes 

using  REBAL  data  for  9  July  1992  1000  CST  unstable  flow  .  19 

2.  Comparison  of  SuperMariah,  Mariah,  and  O’KEYPS  schemes 

using  REBAL  data  for  9  July  1992  2245  CST  stable  flow . 22 


2 


1.  Introduction 


Obukhov  [1]  determined  that  the  sur&ce  boundary  layer  could  be 
characterized  by  a  unique  scaling  length,  L,  whidi  would  represent  the 
surface  layer's  depth.  Monin  and  Obukhov  [2]  generalized  the  Dynamic 
Similarity  of  Flows  Theory  with  the  introduction  of  the  concept  of  the  scaling 
ratio,  z/L,  and  the  log-linear  wind  and  temperature  profiles  for  a  stable 
atmosphere.  Panofsl^  et  al.  [3]  partially  extended  the  similarity  iq)proach 
to  the  unstable  regime.  However,  Panofsky  [4]  provided  the  insight  that  led 
to  the  present  day  linear-quartic  vertical  shapes  of  the  wind,  temperature,  and 
specific  humidity  profiles  in  the  surface  boundary  layer. 

The  linear-quartic  iqpproach  was  dubbed  the  KEYPS*  function  using  the 
initials  of  the  participating  inventors.  Yaglom  [5]  changed  this  acronym  to 
O’KEYPS,  adding  Obukhov's  efforts  of  1946  to  the  original  five  papers. 

Obukhov  did  not  originally  consider  the  effects  of  water  vapor  on 
atmospheric  stability  or  upon  the  structure  of  atmospheric  turbulence.  This 
omission  was  rectified  by  Busch.  [6]  Paulson  [7]  suggested  a  simple 
workable  solution  for  establishing  the  diabatic  influence  functions  but 
Nickerson,  Smiley,  [8]  and  Benoit  [9]  finalized  this  approach. 

Rachele  et  al.  [10]  developed  what  is  known  as  the  Mariah  equations 
whereby  the  Obukhov  [1]  length  could  be  established  from  the  measured 
gradients.  This  simplification  eliminated  the  laborious,  iterative  process  of 
calculating  the  diabatic  influence  functions  to  solve  the  wind,  temperature, 
and  specific  humidity  profile  equations  to  obtain  the  scaling  constants  u*,  6*, 
and  q*  that  are  necessary  components  of  the  surface  energy  balance  equations. 
The  widely  used  aj^roach  is  to  choose  trial  values  of  the  parameters,  then 
cycle  through  an  iterative  process  until  the  solution  converges.  The  intent  of 
this  report  is  to  further  extend  the  Mariah  model  to  what  is  referred  to  as  the 
SuperMariah  approach. 


*Kazanski  and  Monin  (1956),  Ellison  (1957),  Yamamoto  (1959),  Panofoky  (1963),  and  Sellers  (1962). 


3 


2.  Mariah 


The  Mariah  wppro9d\  was  baaed  upon  the  i»emise  that  the  similarity  scaliiig 
constants  for  unstable  flow  can  be  written  as 


u* 


k  Av 
Alnz 


(1) 


e* 


k  Ae 
Alnz 


(2) 


(|>^  A  Inz 


(3) 


where 

u*  =  friction  velocity 

k  -  Karman’s  constant 

V  —  horizontal  wind  speed 

z  =  height  above  surfoce 

0*  =  scaling  temperature 

$  =  mean  potential  temperature 

q*  =  scaling  humidity 

q  =  specific  humidity 

0a  =  cUmensionless  wind  shear 

=  dimensionless  liq)se  rate. 


0a  and  0k  functionally  defined  later.  0^  appears  in  both  equations  Q)  and 
(3)  since  the  similarity  theory  assumes  that  heat  and  mass  transfer  are 
identical. 
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When  evaluating  tfie  above  similarity  constants  we  [11]  found,  using  the 
mean  value  theorem  of  calculus  test,  duU  die  gradients  for  eadi  layer  are 
tangent  to  the  indicated  profile  at  the  height 


z* 


Az 

Afaiz 


instead  of  the  geometric  height  which  is  generally  assumed. 


The  scaling  length  L  is  given  by 


L  = 


(5) 


where 

On  s  virtual  potential  temperature  at  reference  height 

=  6,  (1.0  +  0.61  q)  (6a) 

g  -  acceleration  of  gravity 

6*  —  virtual  scaling  constant  defined  as. 

e;  =  0*  +  0.61  0,  q\  (6t>) 


Substitution  of  equations  (1),  (2),  (3),  (6a),  and  (6b)  into  equation  (S)  yields 

Ow  (Av)* 


L  = 


g  (A0  ^  0.61  01  Lq)  Alnz 


C7) 


as  the  expression  for  L.  For  unstable  (z/L  <  0)  conditions , 


(8a,b) 


and  for  stable  (z/L  >  0)  conditions, 

1  +  P  Y 


(8c) 


Since  in  the  unstable  regime  —  0.^,  equation  (7)  simplifies  to 


L  = 


e,,  (Av)" 

g  (A  6  0.61  Oj  A  q)  A  Inz 


(9) 


In  stable  flow  conditions  where  =  0m  =  (1  +  z/L), 

equation  (7)  becomes 


_ (Av)" _ 

g  0^  (AO  +  0.61  Oj  Aq)  A  Inz 


(10) 


Equation  (10),  after  some  additional  manipulaflon  (see  appendix  A),  becomes 


L  = 


60z*  +  2B 


(11) 


where 

^  2e,,(Av)^ 

g  (A6  +  0.61 6j  Aq)  Alnz 


(12) 


Although  simple  and  efficient,  the  Mariah  parameters  must  be  smoothed  and 
averaged  with  height  in  order  to  invoke  forced  stationarity  upon  the  vertical 
profiles  in  agreement  with  the  tenets  of  similarity.  This  vertical  smoothing 
and  averaging  is  discussed  in  the  next  section. 
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3.  SuperMariah 


SuperMariah  is  an  extension  of  Mariah.  First,  we  determine  if  atmosj^ric 
conditions  are  stable,  unstable,  or  neutral  depending  on  whether  the  potential 
temperature  with  height  is  adiabatic,  superadiabatic,  or  subadiabatic.  After 
determining  the  stability  condition,  we  then  use  the  Mariah  equation  for 
computing  L  using  equation  (9)  for  the  unstable  case  or  equation  (13)  for 
stable  conditions.  Having  L,  we  can  compute  u*,  0*,  and  q*  for  different 
combinations  of  measured  values  of  wind,  temperature,  and  specific  humidity 
from  all  tower  or  mast  levels. 

The  sets  of  u*,  d*,  and  q*  from  all  these  combinations  are  then  averaged 
resulting  in  a  single  estimate  of  each  similarity  constant  for  the  total  profile 
height. 

The  average  values  of  u*,  and  q*  are  inserted  into  equation  (4)  to  obtain 
a  new  value  of  L.  This  estimate  of  L  is  then  used  to  compute  new  values  of 
the  scaling  constants.  This  process  is  repeated  until  all  constants  converge, 
usually  requiring  3  to  8  iterations  based  on  the  following: 

1.  The  Obukhov  length  changes  by  less  than  0.001  m  from  one  iteration  to 
the  next.  In  rare  cases,  the  calculated  value  of  L  decreases  at  each  step  and 
may  become  smaller  than  the  convergence  test,  whereupon  the  results  are 
deemed  invalid. 

2.  An  arbitrary  number  of  steps  has  been  completed  with  no  convergence. 
In  this  case,  results  are  also  invalid. 

Upon  successful  convergence,  the  similarity  profiles  for  wind  speed, 
temperature,  and  specific  humidity  are  calculated  from  the  scaling  constants 
and  Obukhov  length  obtained  (figures  1  and  2).  In  the  following  equations, 
a  subscript  r  indicates  a  value  at  the  reference  height  while  a  subscript  n 
indicates  the  value  at  height  z.. 


For  stable  cases 


(14) 


e 


where 


+  5.0 


L 


For  unstable  cases 


V  =  + 


2  tan'^Cx) 


(15) 


(16) 


(17) 


where 


(18) 


(19) 


(20) 


ai) 
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Figure  2.  SuperMariah  calculated  profiles  compared  to  data  from  REBAL’  92: 
15  July  1»2  2245  CST  stable  flow,  (a)  Temperature  where  *  .0742,  (b)  Wind  speed 
where  u*  =  .184,  (c)  Specific  humidty  where  q*  =  1.474  •  lOe-5. 


4.  The  SuperStar  Option 


E}q)erience  has  shown  that  the  vertical  profile  of  wind,  potential  temperature, 
and  specific  humidity  obtained  by  the  SuperMariah  procure  may  rotated 
slightly  to  the  left  or  the  right  of  the  measured  data  points.  These  profiles 
can  be  adjusted  using  equations  tailored  from  a  paper  by  Rachele  and 
Tunick.  [12]  The  approach  requires  that  adjustment  increments  Av,  Ad,  and 
Aq  be  estimated,  usually  by  taUng  the  difference  between  the  highest  tower 
level  data  values  (or  the  tower  data  level  the  user  has  the  most  confidence  in) 
and  the  corresponding  calculated  value  at  that  level,  z..  These  differences  are 
then  used  in  adjustment  matrices  resulting  in  modifications  to  scaling 
parameters  for  wind,  potential  temperature,  and  specific  humidity,  i.e.,  Au*, 
Ad*,  and  Aq*.  The  procedure  for  one  case  is  demonstrated  in  figure  3. 


The  correction  matrix  for  stable  conditions  is 
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The  correction  matrix  for  unstable  conditions  is 


where 
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Figure  3.  Comparison  of  SupeiMariah  and  SuperStar  profiles,  Nebraska  7  August  lSi56 
1205  CST  (Lettau,  1957).  (a)  Potential  temperature  where  SupeiMariah  #*  =  -.3989, 
Superstar  $*  =  -.4024,  (b)  Wind  spe^  where  SupeiMariah  u*  =  .4741, 
SuperStar  u*  s  .4578,  (c)  Specific  humidity  Miere  SupeiMariah  q*  =  -1.185e-5, 
SuperStar  q*  -  -1.117e>5. 


5.  Discussion 


The  SuperMariah  solution,  altfiou^  iterative,  is  simpler  to  execute  than  the 
O’KEY^  plus  Benoit  [9]  i4)|m>ach  or  the  trial  value  iteration  scheme  using 
the  0*KEYPS  model.  It  also  yields  superior  results,  particularly  if  large 
amounts  of  data  are  being  processed.  Tunick  et  al.  [13]  utilized 
SuperMariah  to  process  the  R^AL  *92  data  sets  (see  Tunick  et  al.  [14]) 
with  outstanding  results. 

The  actual  algorithm  is  quite  simple  and  is  fast  running.  A  comparative 
example  for  unstable  flow  of  Mariah,  SuperMariah  and  the  0*KEYPS 
iq^roach  is  given  in  table  1.  It  must  be  remembered  that  the  0*KEYPS  does 
not  consider  the  effects  of  water  vapor  upon  atmosj^eric  stability  and  is 
based  on  the  premise  that  z/L  »  R  in  an  unstable  atmosirfiere.  llius,  the 
gradient  Richardson  [IS]  numbers  were  first  calculated.  In  order  to  provide 
the  best  estimate  for  the  Monin-C^uldiov  [2]  scaling  ratio,  the  concept  of  the 
i^roximation  of  the  height  derivative  of  R  was  invoked..  Lettau  [16]  was  the 
first  to  suggest  this  approach,  which  may  be  written  as 

t-i  =  ;  I  .  1,2,3,...  (35) 


Both  the  Mariah  and  O’KEYPS  data  were  subjected  to  the  forced  stationarity 
function  of  equation  (35). 

A  similar  comparison  for  the  stable  regime  was  made  and  is  shown  in 
table  2.  For  the  log-linear  0*KEYPS  solution,  it  was  assumed  that 

<|)^  =  1  +  15R,  (36) 


I  15*,) 


(37) 
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with 


(38) 


all  based  iqpon  results  obtained  by  Hansen.  [17] 

The  data  used  to  examine  the  unstable  case  induded  the  REBAL  sanq  ^  for 
KXX)  CST,  9  July  1992.  All  the  calculated  parameten  are  in  agreement  widi 
die  exception  of  die  specific  humidity  scaling  parameter,  q*,  that  shows  up 
as  a  variidile  in  the  Mariah  and  0*KEYPS  solutions.  This  dqmrture  fomi  die 
norm  Q.e.,  a  constant  q*  with  hei^t)  is  due  to  die  REBAL  eiqieriment  being 
conducted  over  bare  soil.  The  evqxvation  rate  over  bare  soil  is  radically 
different  dian  what  would  be  found  over  a  vegetated  surface.  As  a 
consequence,  the  specific  humidity  profiles  and  vertical  gradients  were  not 
entirely  in  equilibrium  with  die  surfece. 

An  examination  of  the  nocturnal  data  of  table  2  for  2245  CST,  15  July  1992, 
indicates  nonequilibiium  specific  humidity  profUes  and  gradients.  The  specific 
hunodity  profile  had  inverted,  a  condition  necessary  for  dew  fall,  but  nddier 
the  relative  humidity  nor  the  gradients  were  sufficient  to  cause  dew. 
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Table  1.  Cbmiiarisoii  of  SuperMariah,  Mariah,  and  0*KEYS  schemes  using  REBAL  for  9  July  1992 
1000  CST  unstable  flow 
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TaUe  1.  Cnnparisoii  ot  SuperMariah,  Mariah,  and  O’KEYS  schemes  using  REBAL  data  tor  9  Jnty  1992 
1000  CST  unstable  flow  (continued) 
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Table  1.  Comparisoii  of  SuperMariah,  Mariah,  and  O’KEYS  schemes  using  REBAL  data  for  9  July  1992 
1000  CST  unstable  flow  (continued) 


Table  2.  Comparison  of  SuperMariah,  Mariah,  and  0*KEYS  srihemcs  using  REBAL  data  fmr  15  July  1992 
2245  CST  unstable  flow  (continueiO 
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6.  Conclusions 


The  methods  presented  in  this  report  are  new,  although  based  on  traditional 
concepts  of  Dynamic  Similarity.  We  have  shown  that  similar  results  may  be 
obtained  by  using  the  O’KEYPS,  Mariah,  or  SuperMariah  i^)proaches.  It  is 
clear  that  the  core  of  the  SuperMariah  methodology  is  based  upon  the 
computation  of  similarity  scaling  constants  using  the  differences  in  the  tower- 
based  measurements  of  wind  speed,  temperature,  and  moisture  (i.e.,  Mariah). 

The  advantages  of  employing  the  SuperMariah  model  are  numerous. 
Primarily,  the  model  executes  quickly.  Secondly,  one  has  the  capability  of 
using  as  many  or  as  few  levels  of  tower  data  as  are  available.  Additionally, 
SuperMariah  outputs  averaged  values  for  the  similarity  scaling  constants 
which  represent  many  subintervals  of  tower  measurements.  Finally,  the 
principle  similarity  premise  of  stationarity  is  preserved.  This  is  too  often 
neglected  in  micrometeorological  tower  data  analyses  schemes. 
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Appendix  A 

Obukhov  Length  for  Stable  Conditions 
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The  Obukhov  lenigth  for  stable  ccmditioiis  can  be  written  as ' 


L 


_ tfw  _ 

g  (A$  *  0.61  Aq)  Alnz 


(A-1) 


or 


gw 

g  (Ad  *  0.61  Aq)  Alnz* 


(A-2) 


From  Hansen* 


(A-3) 


where 


(A-4) 


Substituting  equation  (A-4)  into  (A-3)  and  solving  for  gives 


1  ± 

1+60^ 

id 

1/2 

2 

(A-5) 


*Hinisen,  F.  V.,  1977:  The  Qitical Ridutnbmi Number.  ECOM-5829,  U.S.  Army  Atmoqiheric  Scieaon 
Labotatoty,  White  Sands  Missile  Range,  NM  88002-SS01. 
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Substitutiqg  equation  (A-5)  into  (A-2)  results  in 


L  ±  L(1  ♦  60^)*'* 


2  9^  (Av)» 


g  (Ad  0.61  0^  Aq)  Alnz 
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(A-6) 


or 

±  L(1  +  604^)*'*  =  S  -  L  . 

id 


(A-7) 


Squaring  equation  (A-7)  and  solving  for  L  gives 


L  = 


60z*  +  2B 


(A-8) 


where 


(A-9) 
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